home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / mctst2.f < prev    next >
Text File  |  2002-06-11  |  786b  |  30 lines
  1. C
  2. C PAGE 345-346: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: MCTST2.FOR
  5. C
  6. C NUMERICAL VALUE OF INTEGRAL OVER A DISK IN XY-SPACE BY MONTE CARLO (RANDOM)
  7. C
  8.       DATA  L,N,SUM/256,1000,0.0/     
  9.       F(X,Y) = SIN(SQRT(ALOG(X + Y + 1.0))) 
  10.       PI4 = ATAN(1.0)       
  11.       DO 3 I = 1,N
  12.    2    X = RANDOM(L)       
  13.         Y = RANDOM(L)       
  14.         IF((X - 0.5)**2 + (Y - 0.5)**2 .GT. 0.25)  GO TO 2
  15.         SUM = SUM + F(X,Y)  
  16.         IF(MOD(I,250) .EQ. 0)  THEN   
  17.           VOL = PI4*SUM/REAL(I)       
  18.           PRINT 4,I,VOL     
  19.         END IF    
  20.    3  CONTINUE    
  21.    4  FORMAT(5X,I5,F10.5) 
  22.       STOP
  23.       END 
  24.   
  25.       FUNCTION RANDOM(L)
  26.       L = MOD(16807*L,2147483647)
  27.       RANDOM = REAL(L)*4.6566128752458E-10      
  28.       RETURN
  29.       END 
  30.